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Project  Objectives 


The  project  objective  was  to  significantly  advance  the  quantitative  understanding  of 
solar  and  magnetospheric  energy  inputs  to  the  upper  atmosphere  and  the  resultant 
responses  of  thermospheric  density  and  wind  and  of  ionospheric  electron  density. 
Specifically,  our  goals  were  to: 

•  Quantify  the  distribution  of  magnetospheric  energy  inputs  to  the  thermosphere/ 
ionosphere  system,  and  understand  how  energy  is  partitioned  in  terms  of  generating 
density  and  wind  perturbations  on  local  and  global  scales; 

•  Reveal  and  understand  new  aspects  of  the  thermosphere  response  to  solar  extreme 
ultraviolet  (EUV)  variability,  and  how  this  response  interacts  with  that  to  magnetospheric 
inputs; 

•  Develop  improved  parameters  for  first-principles  models  of  the 
thermosphere/ionosphere  system,  in  order  to  optimize  their  ability  to  simulate  the 
system’s  response  to  variable  energy  inputs. 

Work  Carried  Out  and  Results  Obtained 

1 .  Seasonal  variation  of  thermospheric  density  and  composition. 

Qian  et  al.  [2008]  investigated  the  causes  of  seasonal  variations  in  thermospheric 
neutral  density  and  composition,  which  exhibit  maximum  densities  near  the  equinoxes, 
a  primary  minimum  during  northern  hemisphere  summer,  and  a  secondary  minimum 
during  southern  hemisphere  summer.  This  pattern  of  variation  is  described  by 
thermospheric  empirical  models.  However,  the  mechanisms  are  not  well  understood. 
The  annual  insolation  variation  due  to  the  Sun-Earth  distance  can  cause  an  annual 
variation;  large-scale  inter-hemispheric  circulation  can  cause  a  global  semiannual 
variation;  and  geomagnetic  activity  can  also  have  a  small  contribution  to  the  semiannual 
amplitude.  However,  simulations  by  the  NCAR  Thermosphere-Ionosphere- 
Electrodynamics  General  Circulation  Model  (TIE-GCM)  indicate  that  these  seasonal 
effects  do  not  fully  account  for  the  observed  annual/semiannual  amplitude,  primarily  due 
to  the  lack  of  a  minimum  during  northern  hemisphere  summer.  A  candidate  for  causing 
this  variation  is  change  in  composition,  driven  by  eddy  mixing  in  the  mesopause  region. 
Other  observations  and  model  studies  suggest  that  eddy  diffusion  in  the  mesopause 
region  has  a  strong  seasonal  variation,  with  eddy  diffusion  larger  during  solstices  than 
equinoxes,  and  stronger  turbulence  in  summer  than  in  winter.  Qian  et  al.  determined  a 
seasonal  variation  of  eddy  diffusion  that  is  compatible  with  this  description.  Simulations 
showed  that  when  this  function  is  imposed  at  the  lower  boundary  of  the  TIEGCM, 
neutral  density  variation  consistent  with  satellite  drag  data,  and  O/N2  consistent  with 
measurements  by  TIMED/GUVI,  are  obtained.  These  model-data  comparisons  and 
analyses  indicate  that  transmission  of  turbulent  mixing  from  the  lower  atmosphere  may 
contribute  to  seasonal  variation  in  the  thermosphere,  particularly  the  asymmetry 
between  solstices  that  cannot  be  explained  by  other  mechanisms. 


2.  Empirical  model  of  high-latitude  electrodynamics. 

We  developed  empirical  models  of  high-latitude  electric  potential,  magnetic  stream 
function,  electric-field  variability,  and  downward  Poynting  flux,  based  on  measurements 
from  the  Dynamics-Explorer-2  spacecraft.  Each  model  is  a  function  of  magnetic  latitude 
and  magnetic  local  time,  and  is  parameterized  in  terms  of  hemisphere,  dipole  tilt  angle 
toward  the  Sun  (season),  and  magnitude  and  orientation  of  the  interplanetary  magnetic 
field  (IMF).  The  models  of  electric  potential  and  magnetic  stream  function  are  similar  to 
comparable  models  by  Weimer  [Weimer,  D.R.  (2005),  Improved  ionospheric 
electrodynamic  models  and  application  to  calculating  Joule  heating  rates,  J.  Geophys. 
Res.,  110,  A05306,  doi:10.1029/2004JA010884],  but  the  models  of  electric-field 
variability  and  Poynting  flux  contain  important  additional  information  related  to  the 
variable  components  of  the  electric  and  magnetic  fields  about  their  climatological 
values.  In  particular,  the  variable  field  components  produce  an  important  additional 
energy  input  in  the  day-side  high-latitude  region.  The  empirical  models  are  designed  for 
use  in  simulation  models  of  upper-atmospheric  dynamics,  as  representations  of  the 
magnetospheric  forcing  and  electromagnetic  energy  input.  They  have  been  used  by 
Deng  et  al.  [2009]  to  drive  the  TIE-GCM  (see  task  3).  The  component  models  have 
been  prepared  as  computer  subroutines  that  will  be  made  available  to  the  community 
after  completion  of  testing  and  documentation.  Examples  of  the  fields  predicted  by 
these  models  over  the  northern  polar  region  are  shown  in  Figures  1-12  at  the  end  of  this 
report. 

3.  Impact  of  electric-field  variability  on  the  thermosphere. 

Deng  et  al.  [2009]  used  the  empirical  models  of  electric  potential,  electric-field 
variability,  and  downward  Poynting  flux  described  in  task  2  as  magnetospheric  driving 
inputs  to  the  TIE-GCM,  to  test  how  well  the  total  Joule  heating  produced  by  the 
combined  mean  and  variable  components  of  the  electric  field  agree  with  the 
independent  estimate  of  this  heating  given  by  the  Poynting  flux.  Whereas  the  TIE-GCM 
values  of  Joule  heating  depend  on  the  ionospheric  conductivity,  the  Poynting  flux  does 
not,  and  so  these  are  two  independent  measures  of  the  energy  input  to  the 
thermosphere.  Although  some  differences  in  spatial  patterns  of  the  Joule  heating  and 
Poynting  flux  were  noted,  the  hemispherically  integrated  values  agreed  quite  well, 
demonstrating  that  inclusion  of  the  electric-field  variability  significantly  improves  our 
ability  to  estimate  magnetospheric  energy  inputs  to  the  thermosphere. 

4.  Dependence  of  the  high-latitude  thermospheric  densities  on  the  interplanetary 
magnetic  field. 

Kwak  et  al.  [2009]  carried  out  a  systematic  analysis  of  the  interplanetary  magnetic  field 
(IMF)  By  and  Bz  influences  on  observed  thermospheric  density.  For  this  purpose,  the 
high-latitude  southern  summer  thermospheric  total  mass  density  near  400  km  altitude, 
derived  from  the  high-accuracy  accelerometer  on  board  the  Challenging  Minisatellite 
Payload  (CHAMP)  spacecraft,  was  statistically  analyzed  in  magnetic  coordinates.  The 
difference  density  distributions,  obtained  by  subtracting  values  for  zero  IMF  from  these 


for  nonzero  IMF,  vary  strongly  with  respect  to  the  direction  of  the  IMF:  Difference 
densities  for  negative  By  show  significant  enhancements  in  the  early  morning  hours  and 
hours  around  dawn,  but  show  reduced  values  in  the  dusk  sector.  For  positive  By  the 
difference  densities  are  opposite  in  sign  to  those  for  negative  By.  Density  differences  for 
negative  Bz  show  significant  increases  in  the  cusp  region  and  premidnight  sector,  but  a 
small  decrease  in  the  dawn  sector.  The  difference  densities  at  high  latitudes  tend  to  be 
weakest  when  Bz  is  positive.  The  authors  suggested  that  the  high-latitude  thermospheric 
density  variations  for  different  IMF  conditions,  especially  in  the  dawn  and  dusk  sectors, 
can  be  strongly  determined  by  thermospheric  winds,  which  are  associated  with  the 
ionospheric  convection  and  vary  strongly  with  respect  to  the  IMF  direction.  They  also 
suggested  that  density  variations,  especially  in  auroral  and  cusp  regions,  are  also 
influenced  by  the  local  heating  associated  with  ionospheric  currents,  which  vary  with 
IMF  conditions. 

5.  Thermospheric  density  enhancements  in  the  dayside  cusp  region  during  strong  IMF 
By  conditions. 

Crowley  et  al.  [2010]  provided  the  first  detailed  explanation  of  a  high-latitude  density 
enhancement  in  the  cusp  sector  observed  in  tri-axial  accelerometer  data  from  the 
Challenging  Minisatellite  Payload  (CHAMP)  satellite,  focusing  on  the  August  24,  2005 
interval.  The  Thermosphere  Ionosphere  Mesosphere  Electrodynamics  General 
Circulation  Model  (TIMEGCM)  was  driven  by  high-fidelity  high-latitude  inputs  specified 
by  the  Assimilative  Mapping  of  Ionospheric  Electrodynamics  (AMIE)  algorithm,  and 
reproduced  the  main  features  of  the  density  enhancements.  The  TIMEGCM  and  AMIE 
provide  a  global  framework  for  interpretation  of  the  CHAMP  densities.  The  simulations 
reveal  that  the  observed  density  enhancement  in  the  dayside  cusp  region  resulted  from 
unexpectedly  large  amounts  of  energy  entering  the  ionosphere-thermosphere  system  at 
cusp  latitudes  during  an  interval  of  strong  (+20  nT)  By. 

6.  A  computationally  compact  representation  of  Magnetic-Apex  and  Quasi-Dipole 
coordinates  with  smooth  base  vectors. 

Many  structural  and  dynamical  features  of  the  ionized  and  neutral  upper  atmosphere 
are  strongly  organized  by  the  geomagnetic  field,  and  are  most  conveniently  represented 
with  respect  to  magnetic  coordinates.  Modified  Apex  coordinates  are  appropriate  for 
calculations  involving  electric  fields  and  magnetic-field-aligned  currents;  Quasi-Dipole 
coordinates  are  appropriate  for  calculations  involving  height-integrated  currents. 

Emmert  et  al.  [2010]  developed  a  compact,  smooth,  and  robust  representation  of  the 
transformation  from  geodetic  to  Quasi-Dipole  and  Modified  Apex  coordinates.  They 
made  the  computer  code  available  to  the  community  for  evaluating  the  coordinates  and 
base  vectors,  as  auxiliary  material  to  their  publication. 

7.  Model  simulation  of  thermospheric  response  to  recurrent  geomagnetic  forcing. 

Qian  et  al.  [2010]  assessed  model  capability  in  simulating  thermospheric  response  to 
recurrent  geomagnetic  forcing  driven  by  modulations  in  the  solar  wind  speed  and  the 
interplanetary  magnetic  field.  Neutral  density  and  nitric  oxide  (NO)  cooling  rates  were 


simulated  for  the  declining  phase  of  solar  cycle  23.  The  simulated  results  are  compared 
to  neutral  density  derived  from  satellite  drag  and  to  NO  cooling  measured  by  the 
Thermosphere  Ionosphere  Mesosphere  Energetics  and  Dynamics  (TIMED)  sounding  of 
the  atmosphere  using  broadband  emission  radiometry  (SABER)  instrument.  Model-data 
comparisons  show  good  agreement  between  the  model  and  the  measurements  for 
multiday  oscillations,  as  well  as  good  agreement  for  longer-term  variations.  The 
simulations  demonstrate  that  the  multiday  oscillation  of  density  is  globally  distributed  in 
the  upper  thermosphere  but  restricted  to  high  latitudes  in  the  lower  thermosphere.  The 
density  variation  in  the  upper  thermosphere  exhibits  less  latitude  dependence  than  the 
temperature  variation  because  of  the  effects  of  composition  changes.  Model  simulations 
also  show  that  NO  density  and  temperature  play  primary  roles  in  the  multiday  oscillation 
of  NO  cooling  rates. 

8.  On  the  ionospheric  application  of  Poynting’s  Theorem. 

It  has  been  proposed  that  the  geomagnetic-field-aligned  component  of  the  perturbation 
Poynting  vector  above  the  ionosphere,  as  obtained  from  the  cross  product  of  the  electric 
and  magnetic-perturbation  fields  observed  on  a  spacecraft,  may  be  used  to  estimate  the 
field-line-integrated  electromagnetic  energy  dissipation  in  the  ionosphere  below. 
Richmond  [2010]  clarified  conditions  under  which  this  approximation  may  be  either  valid 
or  invalid.  It  was  shown  that  the  downward  field-aligned  component  of  the  perturbation 
Poynting  vector  can  underestimate  the  electromagnetic-energy  dissipation  in  regions  of 
high  ionospheric  Pedersen  conductance,  and  it  can  significantly  overestimate  the 
dissipation  in  regions  of  low  conductance.  Local  values  of  upward  perturbation  Poynting 
vector  do  not  necessarily  correspond  to  net  ionospheric  generation  of  electromagnetic 
energy  along  that  geomagnetic-field  line.  An  Equipotential-Boundary  Poynting-Flux 
(EBPF)  theorem  was  presented  for  quasi-static  electromagnetic  fields  as  follows:  when 
a  volume  of  the  ionosphere  is  bounded  on  the  sides  by  an  equipotential  surface  and  on 
the  bottom  by  the  base  of  the  conducting  ionosphere,  then  the  area  integral  of  the 
downward  normal  component  of  the  perturbation  Poynting  vector  over  the  top  of  that 
volume  equals  the  energy  dissipation  within  the  volume.  This  equality  does  not  apply  to 
volumes  with  arbitrary  side  boundaries.  However,  the  EBPF  Theorem  can  be  applied 
separately  to  different  components  of  the  electric  potential,  such  as  the  large-  and 
small-scale  components.  Since  contours  of  the  small-scale  component  of  potential  tend 
to  close  over  relatively  localized  regions,  the  associated  small-scale  structures  of 
downward  perturbation  Poynting  vector  tend  to  be  dissipated  locally. 

9.  Electrodynamics  of  ionosphere-thermosphere  coupling. 

Richmond  [2011]  presented  an  overview  of  ionosphere-thermosphere  electrodynamic 
coupling.  Collisions  between  the  charged  and  neutral  constituents  of  the  upper 
atmosphere  couple  their  respective  dynamics  and  energetics.  Magnetic  stresses  readily 
transfer  momentum  and  energy  over  long  distances  along  geomagnetic-field  lines, 
accompanied  by  electric  fields  and  currents.  Consequently,  the  E  and  F  regions  of  the 
ionosphere  are  strongly  coupled,  and  momentum  is  transferred  between  the  lower  and 
upper  thermosphere  through  the  currents  and  their  associated  ion  drag.  Electrical 


conductivity  mediates  the  degree  of  ion-neutral  coupling.  Conductivity  is  highly  variable, 
and  is  itself  affected  by  the  electric  field  in  various  ways.  Thermospheric  winds  drive  the 
ionospheric  wind  dynamo.  The  winds  are  created  by  daily  absorption  of  solar  radiation 
in  the  thermosphere,  by  upward  propagating  solar  and  lunar  tides,  by  ion-drag 
acceleration  at  high  latitudes,  and  by  Joule  heating  at  high  latitudes.  Electric  current 
flows  globally  in  the  ionosphere  and  along  geomagnetic-field  lines  through  the 
magnetosphere.  Interactions  between  the  ion  and  neutral  motions  produce  feedbacks 
that  affect  the  dynamics  of  both  components.  Simulation  models  of  thermosphere- 
ionosphere-electrodynamic  interactions  provide  powerful  tools  for  investigating  the 
nature  of  these  interactions,  and  for  testing  how  well  the  uncertain  model  inputs  and  the 
physics  incorporated  in  the  models  are  able  to  predict  observed  features  of  the 
ionosphere  and  thermosphere. 

10.  Thermospheric  mass  density:  an  overview  of  temporal  and  spatial  variations. 

Qian  and  Solomon  [2011]  summarized  and  discussed  thermospheric  density  variations, 
their  magnitudes,  and  their  forcing  mechanisms,  using  neutral  density  data  sets  and 
modeling  results.  The  density  data  sets  included  observations  by  accelerometers  on  the 
Challenging  Mini-satellite  Payload  (CHAMP)  and  measurements  of  orbit  changes  on 
other  space  objects  due  to  drag  at  perigee,  and  global-mean  density  derived  from 
thousands  of  orbiting  objects.  Modeling  results  are  from  the  TIE-GCM  and  from  the 
NRLMSISE-00  empirical  model.  Neutral  density  shows  complicated  temporal  and 
spatial  variations  driven  by  external  forcing  of  the  thermosphere/ionosphere  system, 
internal  dynamics,  and  thermosphere  and  ionosphere  coupling.  Temporal  variations 
include  abrupt  changes  with  a  time  scale  of  minutes  to  hours,  diurnal  variation,  multi¬ 
day  variation,  solar-rotational  variation,  annual/semiannual  variation,  solar-cycle 
variation,  and  long-term  trends  with  a  time  scale  of  decades.  Spatial  variations  include 
latitudinal  and  longitudinal  variations,  as  well  as  variation  with  altitude.  Atmospheric 
drag  on  satellites  varies  strongly  as  a  function  of  thermospheric  mass  density.  Errors  in 
estimating  density  cause  orbit  prediction  error,  and  impact  satellite  operations,  including 
accurate  catalog  maintenance,  collision  avoidance  for  manned  and  unmanned  space 
flight,  and  re-entry  prediction. 
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Figure  1 .  Electric  potential  as  a  function  of  IMF  clock  angle  for  a  transverse  IMF 
magnitude  of  5  nT  and  for  zero  dipole  tilt. 
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Figure  3.  Electric  potential 


as 
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function  of  dipole  tilt  for  a 
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southward  IMF  of  5  nT. 
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Figure  4.  Standard  deviation  of  the  northward  electric  field  as  a  function  of  IMF  clock 
angle  for  a  transverse  IMF  magnitude  of  5  nT  and  for  zero  dipole  tilt. 
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Figure  5.  Standard  deviation  of  the  northward  electric  field  as  a  function  of  southward 
IMF  magnitude  for  zero  dipole  tilt. 
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Figure  6.  Standard  deviation  of  the  northward  electric  field  as  a  function  of  dipole  tilt  for 
a  southward  IMF  of  5  nT. 


MagPot 

Btrans-  5  0  rT 
Dipole  tilt  angle-  0* 


Figure  7.  Magnetic  stream  function  as  a  function  of  IMF  clock  angle  for  a  transverse  IMF 
magnitude  of  5  nT  and  for  zero  dipole  tilt. 
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Figure  8.  Magnetic  stream  function  as  a  function 
dipole  tilt. 
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Figure  9.  Magnetic  stream  function  as  a  function  of  dipole  tilt  for  a  southward  IMF  of  5 


nT. 
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Figure  10.  Downward  Poynting  flux  as  a  function  of  IMF  clock  angle  for  a  transverse 
IMF  magnitude  of  5  nT  and  for  zero  dipole  tilt. 
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Figure  1 1 .  Downward  Poynting 
dipole  tilt. 
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Figure  12.  Downward  Poynting  flux  as  a  function  of  dipole  tilt  for  a  southward  IMF  of  5 
nT.  (The  top  plot  is  for  -24°  tilt,  the  lower-left  plot  is  for  0°,  and  the  bottom-right  plot  is  for 
+24°  tilt.) 


